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Quantitative entanglement witnesses allow one to bound the entanglement present in a system 
by acquiring a single expectation value. In this paper we analyze a special class of such observ- 
ables which are associated with (generalized) Werner and Isotropic states. For them the optimal 
bounding functions can be easily derived by exploiting known results on twirling transformations. 
By focusing on an explicit local decomposition for these observables we then show how simple clas- 
sical post-processing of the measured data can tighten the entanglement bounds. Quantum optics 
implementations based on hyper-entanglement generation schemes are analyzed. 
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I. INTRODUCTION 

Detecting and characterizing entanglement in quantum 
systems is a major issue of quantum information the- 
ory, both from the theoretical and from the experimental 
point of view [IIH]. In this context Entanglement Wit- 
nesses (EWs) turn out to be useful tools that allow one to 
address this problem with a minimal experimental effort, 
even in those cases where a complete tomographic recon- 
struction of the state is not experimentally accessible [3] . 
Indeed, this technique merely requires one to measure 
a single (non local) observable whose expectation value 
yields a certifiable evidence of the presence of entangle- 
ment in the system. Unfortunately it cannot provide a 
definitive answer for all system configurations (indeed in 
most of the cases, the result of the measurement will be 
inconclusive). Furthermore, even when successful, the 
result of an EW measure will not be useful in general 
to quantify the amount of entanglement present in the 
system. Still it was recently pointed out that the 

outcome of such measurements can be used to provide 
at least non trivial (lower) bounds for the latter. These 
observations led to the notion of Quantitative Entangle- 
ment Witness (QEW) [S], which identify those operators 
whose expectation values can be translated into a quan- 
tification of the entanglement present in the measured 
system. 

From a purely technical point of view, it is worth stress- 
ing that even though QEWs were originally discussed in 
the context of standard EWs, they can be thought as 
independent from the latter. Indeed an EW is character- 
ized by well established properties, e.g. it must assume 
positive values on all separable states, and it must have 
at least one negative eigenvalue On the contrary 
the definition of QEW can be applied to any operator Q 
for which it is possible to associate a non trivial hound- 
ing function (see next section) for a given entanglement 
measure. Notable examples of such QEWs were already 
known in literature even before the formal introduction 
of this notion, e.g. see [Hj. In particular the idea of 



bounding the entanglement of a state via the acquisition 
of a single expectation value was addressed in a series of 
papers [1^ that applied previous ideas on twirling trans- 
formations [TT] and Werner states [HI [13]; the same argu- 
ment has been extended lately to Isotropic states [HI [T51 - 
[T9] . The geometrical properties of Werner and Isotropic 
states can be used to build explicitely QEW capable of 
detecting up to a high degree of bipartite entanglement. 
In this paper we review the derivation of such QEW |S]-[S] 
adopting a formal yet simple scheme. Subsequently us- 
ing an explicit decomposition of such operators in terms 
of correlated local observables [20] we perform a data- 
processing optimization which will enable further tight- 
ening of the entanglement bounds. Expressing EW or 
QEW in terms of local operation and classical communi- 
cation (LOCC) measurement schemes is yet another im- 
portant aspect of entanglement detection [1, 9, 19 26 . It 
deals with the problem of acquiring expectations values 
of a non-local observable (i.e. an observable that does not 
admit a separable basis of eigenvectors) , when only local 
resources are accessible experimentally. Even though the 
decomposition we assume is not the most efficient one in 
terms of setup preparations [T], it has the main advan- 
tage of being easily implementable in some specific exper- 
imental configurations which exploit multi-rail quantum 
encoding strategies [?7lI33j . We will explicitly focus on 
some optical scheme based on hyper- entanglement setups 
[?7r.29J where bipartite entangled qudits states are gen- 
erated by exploiting spatial and polarization degrees of 
freedom. 

The paper is organized as follows. After reviewing 
some basic definitions around quantitative entanglement 
witness, in section |lT] we focus on QEWs of Werner class 
and of Isotropic class for which we present the optimal 
bounding functions. In Sec. jlllj we introduce the local 
decompositions for such QEWs which use 2d indepen- 
dent set-up preparations (d being the local dimension of 
a system site). We then show how the same measure- 
ment outcomes can be post-processed to strengthen the 
entanglement bounds. Sec. |IV| finally discusses how to 
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implement the proposed measurement scheme in the con- 
text of optical multi-rail encodings. 

II. QUANTITATIVE ENTANGLEMENT 
WITNESSES 

Let Q be an observable acting on a bipartite system 
AB, and let £m be a given convex and continuous mea- 
sure of bipartite entanglement, e.g. the entanglement of 
formation [H El [U]. We say that Q is a QEW for the 
measure £m if can define a (nontrivial) non-negative 
function b for which it holds 



bound on £m (p) j namely [1] 

SMip)> max 6opt (Tr [T* (Q) p] ) = 6opt {q* ) , (6) 

where the maximization is now performed over the set of 
all LOCC unitaries, and where g* indicates the value of 
Tr [r* (Q) p] that yields the highest value of 6opt . Unfor- 
tunately q.^, is a rather complicated (nonlinear) function 
of the state p which is arguably hard to compute. Still, 
we will see in the following that, for the cases we con- 
sider here, an operational characterization of it can be 
obtained by focusing on a subclass of F. 



Sm (p) > b (Tr 



(1) 



A. Werner and Isotropic QEWs 



for all p density matrices of AB. Analogously to a stan- 
dard entanglement witness [IHl] , a QEW provides an op- 
erational method for detecting entanglement. Further- 
more it also yields a quantitative estimation of it. As 
pointed out in [SJ O [H] j ^ fundamental problem is to find 
the optimal bounding function of a given QEW Q, i.e. 
the highest positive b for which ([T]) holds. The latter can 
be formally defined as 



"opt I 



ix)^mi{£Mip)\Tr[Qp]^x} 



(2) 



where the minimization is performed over the whole set of 
density matrices. This is a convex (generally not mono- 
tonic) function which typically is rather difficult to ob- 
tain. A partial characterization of 6opt was proposed in 
[6] where it was pointed out that, due to the convexity 
of entanglement measures, we can write 



bopt{x) = sup {Xx - Cm 

A 



(3) 



where the supremum is taken over all real values of the 
parameter A and where Cm is the Legendre transform 
of the measure functional £m i with respect to the trace 
scalar product, i.e. 



Cm {A) = sup {Tr [A p] - £m {p)) 
p 



(4) 



(here the supremum is taken over the whole set of density 
matrices). 

A useful property of the set of QEW is that it is invari- 
ant under Local Operation and Classical Communication 
(LOCC). In particular any (non trivial) unitary LOCC 
transformation F takes a given QEW into a new QEW 
while preserving its bounding functions including the op- 
timal one, which is trivially again optimal |34j . This is 
a simple consequence of the cyclicity of the trace which 
gives 



£M[p)^£M{T{p))>b{Tr[T*m p]).. 



(5) 



for all 6 of Q (in this expression F* stands for the adjoint 
super-operator of F). Notice that exploiting the proper- 
ties of the transformations F, we can derive a stronger 



Finding the optimal bounding functions for QEWs is a 
hard problem [51 El ^ . However there are certain highly 
symmetric choices of Q for which 6opt can be obtained 
without going through the optimizations of (|3|. Proto- 
typical examples are given by the operators 



F= ^ |a/3)(/3c 



= ^ |aa)(/?/3|, (7) 



a, 13=1 



a, 13=1 



where {|Q;);a — \,---d\ is the canonical basis for the 
subsystem A and B (both assumed to be d-dimensional) . 
The first is nothing but the swap operator on the compos- 
ite system AB, while the second is a projector into the 
(not normalized) maximally entangled state with con- 
stant coefficients. They are related via partial transpose 
(i.e. F = G'^) and are QEWs which we can call of Werner 
class and of Isotropic class respectively, for reasons which 
will be clear soon. Their optimal bounding functions can 
be computed for several entanglement measures includ- 
ing entanglement of formation {£of)i relative entropy of 
entanglement |13) . and convex-roof extended negativity 
|18j . For the sake of simplicity however in the following 
we focus only on the former. In this case the optimal 
bounding functions for F and G are respectively given by 



OWcr 



{x) 




blso{x) = 



co[/i2(7) 




for X < 0, 
for X > 0, 



(l-7)log(d- 



(8) 



1)] for a; > 1, 
for a; < 1, 



where h2{y) = —ylogy — (1 — y)log(l — y) is the bi- 
nary Shannon entropy, 7 is the function 7(0;) = -^{y/x + 
\/ {d — l){d ~ x)), and 'co' means taking the convex- hull 
of the inner expression, i.e. the largest convex curve 
nowhere larger than the given one - notice however that 
in Q the argument is already convex when d = 2 [35 . 

Proving that the operator F is a QEW with optimal 
bounding function ([s]) involves the knowledge of the prop- 
erties of Werner states [T^J [T3] and notion of twirling 
transformations [TT] . Werner states pw in an T^d ^ Hd 
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space are those which are unvariant under the whole 
U (E)U group of transformations 



{U (E)U)pw{U (E)U)^ = pw, 



(9) 



for all unitaries U on Jid- They form a one-parameter 
manifold of states which can expressed as 



Pwif) 



1 



d(d2 - 1 



[(d-/)l + (/d-l)F], (10) 



where 1 is the (P x cP identity operator and where F 
is the swap operator ([t]). The parameter / takes values 
between -1 and 1 to ensure positivity, and it coincides 
with the expectation value of F over pw[f ), i-e. 



/ = Tr[Fpw(/)] 



(11) 



The entanglement of formation Eqf {pwif)) for state 
Pwif) is a known quantity which can be expressed 
in terms of the function (|8|, i.e. 

(pwif)) - fewcr(/) = ^'Wcr (Tr [¥ pw]) ■ (12) 

To demonstrate that bwor is a valid bounding function 
for F let us define the following completely-positive trace- 
preserving mapping, acting on a generic density matrix 
p of AB, 

Va{p)= I dp{U){U®U)p{U®U)\ (13) 



where A is the unitary group of "Hd and dp,, its Haar 
measure, is uniquely defined by compactness of A. The 
transformation T'a is often labeled as twirling operation 
[llj . Three properties of Pa are important to mention 
for our purposes [HI [13l [18] : first, Vk is a projector, i.e. 

= Va, and its range is exactly the ID manifold of 
Werner states. Moreover, Va preserves the expectation 
value of F, i.e. Tr [¥ p] = Tr [F7'a(p)], which, along with 
( 11 1, tells us that 



rA{p)^Pw (Tr[Fp]) 



(14) 



Finally, since its Kraus decomposition ( 13 ) is made of 



tensor product operators, Va is a (non-invertible) LOCC 
transformation. Therefore it can only decrease entangle- 
ment, thus providing the following inequality 



^f(Pa(p)) < SoFip). 



(15) 



Putting together (12), ( (m] ) and (15) wc finally obtain 

[iniiniiiiiiH] 



53f(p)> fewer (Tr[Fp]), 



(16) 



for any given state p, which is the proper definition of 
bounding function for fewcr- To check optimality, it is 
sufficient to see that for any x E [—1,1] there is a density 
matrix p with Tr [¥ p] = x for which equality holds: this 
happens for the Werner state pwix). Then, since for any 



function b higher in x than fcwcr, the condition ([!]) would 
be broken, 6wcr must be optimal. 

The proof of the optimality of the function ^ for the 
Isotropic QEW G is similar to the one we exploited in 
the Werner case. This time we consider isotropic states 
P^HTT] . i.e. the set of states p/ invariant under every 
transformation of the form 



{u®u*)piiu®u*y = pi, 



(17) 



where "*" denotes complex conjugation of the matrix en- 
tries, say, with respect to the canonical basis. Isotropic 
states form again a one-parameter family, defined by 



Piig) 



1 



did^ - 1) 



[{d-g)l + igd-l) 



(18) 



with < g < d and G as in ([7). The parameter g is also 
the expectation value g — Tr Grpi{g)]. On top of that, 
their entanglement of formation is known |15| to be given 
by Sopipiig)) — bisoig)- Again, we exploit the properties 
of a twirling operation, this time of the form 



V'a{p)= / dp{U) {U ®U*)p{U ®U*)^ 



(19) 



which projects [131 118) any given state p in the isotropic 
state p/(Tr [G p\) while preserving the expectation value 
of G and degrading the entanglement of formation. Fol- 
lowing these considerations we get [TH] 



£of{p) > foho(Tr[Gp]), 



(20) 



so that 6iso is a bounding function for G, and it is opti- 
mal since any function pointwise higher would break the 
bounding condition ([T]) upon one isotropic state at least 
(it is worth pointing out that the inequality ( 20 ) was 



originally derived in [9 without invoking the properties 
of twirling). 



III. LOCAL DECOMPOSITION OF QEW 

A practical problem in entanglement detection is how 
to acquire expectation values of non-local observables (for 
instance our QEWs), while we are typically allowed to 
perform only LOCC measurements [TJ EOH^S] . The stan- 
dard solution is to expand the witness operator Q into a 
sum of local observables of the form 



B, 



(21) 



with Ci being some complex coefficients. The expecta- 
tion value of Q can then be reconstructed by simple data 
post-processing from the expectation values of the ten- 
sor product observables Ai (g) Bi which can be acquired 
by performing classical correlated measures on A and B. 
The decomposition (21) is not unique and, unless Q ad- 
mits a set of separable eigenvectors [Mj, it will involve a 
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certain number of non-commuting operators Aii^ Bi. In 
this case several independent experimental setups need 
to be prepared, one for each non-commuting set of ob- 
servables Ai (g) Bi entering into (21). It is thus reason- 



able to try and optimize the local decomposition in a 
way that it corresponds to the smallest possible experi- 
mental effort. This problem gathered a lot of attention 
POI [?n [21H1S] but it appears that no general solution 
has been found yet. Instead, several specific solutions for 
particular classes of witnesses have been proposed. 

In this section we focus on a local decomposition of the 
QEWs of the Werner and Isotropic class which was first 
discussed in [50]. Such solution is not the most efficient 
one in terms of number of setup re-preparation. Indeed 
it requires ~ 2d independent setups, while a theoretical 
lower bound of d was presented in [26, which, at least for 
d is prime is known to be achievable [25j . Yet the pro- 
posed decomposition presents at least two advantages. 
First, as we shall see in Section [lII A[ by simply reprocess- 
ing the data accumulated during the various experiment 
runs, one can obtain expectation values of an (exponen- 
tially large) class of QEW operators which are LOCC- 
equivalent in the sense of ^ to the one we started with. 
This leads us to the experimental acquisition of several 
entanglement bounds which are optimal in the sense of 
([2|: upon taking their highest value we can hence pro- 
vide an estimation of the quantity of ([6]). Secondly, 



as we shall see in Section IV the proposed decomposi- 
tion scheme appears to be naturally suited for multi-rail 
encoding implementations. 

To start we introduce the following set of (local) single- 
site observables 



Pa \a){a\ , 
Yap = i\a}{f3\-i\l3){a\, 



(22) 



defined for any a, /3 S {0, • • • , d}, a ^ /3 (as before the 
vectors \a) stands for the canonical basis of the site). 
Since they form a basis for the space of operators acting 
on 7id, we can express F and G as follows. 



F 



G 



a 



Pc^ + jYl ^^°^f> ® + Yap ® Yaf}) ., 



-E 



{Xap «) Xap - Yap (E) Y^p) ■ 



^ (23) 
which is of the form (21) with terms and with the 



coefficients being real. Therefore the expectations val- 
ues of F and G can both be acquired with at most order 
d^ independent experimental setups. Of course, some of 
the operators appearing in the sums (23) commute with 



each other, so in principle it is possible to evaluate them 
within the same experiment. For instance any pair of op- 
erators of the set S = {Pa, Xap, Yap} commute iff all 
the indexes are different, and clearly this commutation 
rule extends to all the observable of the form A^A with 





Figure 1: (colour online) Sketch of the minimization strategy 
for d = 7, first two steps. A red line connecting vertices 
a and /3 means measuring Xa/i ® Xa/s, the green dot is the 
measurement of <8> Pa • 



A e S. Exploiting this simple property one can indeed 
reduce |20) the number of setup re-preparations needed 
to obtain the all values {Pa <S) Pa), {Xap <E> Xap) and 
{Xap ® Xap) from d^ to order d (here we use (• • • ) to 
represent the expectation value with respect to the state 
p). For the sake of clarity, we review this result by re- 
deriving using a different method than that presented in 
[20] . To do so we find it convenient to turn the prob- 
lem into a graph colouring game. Let us assume that 
we have a graph made of d vertices and all the possible 
edges among them (also known as complete graph), e.g. 
see figures [l] and [2] we are also given one green brush to 
colour the vertices with, alongside with one blue and one 
red brushes to paint the edges of the graph. The idea is 
that when we paint the a-th vertex green we are actu- 
ally acquiring {Pa ® Pa), while if we paint in red (resp. 
blue) the edge connecting the nodes a and /3 we evaluate 
{Xap ® Xap) (resp. {Yap ® Yap))- As a rule, during a 
single turn the player can paint everything he wants, but 
not: 

• an edge and a vertex at one of its ends, 

• two consecutive edges, 

• the same edge twice; 

clearly, these game rules encode the commutation rela- 
tions we discussed. The aim of this game is to paint all 
the vertices green and all the segments in both red and 
blue spending the minimal number of turns. Optimal 
solution for the game are presented next: 

Colouring for d odd - For simplicity, we draw the graph 
like a regular polygon with d edges, including all the di- 
agonals, see figure [T] At the first turn the player paints 
in red one edge and all the diagonals that are parallel to 
that edge, as well as the farthest vertex, obviously green. 
From turn 2 to d he chooses everytime another edge of 
the polygon and repeats the process, until all the edges 
have been chosen. From turn d+1 to 2d he does the same 
thing he did in the first d turns but using the blue brush 
instead of the red one. That completes the game, in 2d 
turns. This colouring strategy is optimal because, due 
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Figure 2: (colour online) Sketch of the minimization strategy 
for d = 6. Turns 1, 2 and 11 are drawn. The six-th vertex 
has been placed in the geometrical center of the pentagon. 



to the rules, every turn a player cannot paint more than 
{d — l)/2 segments, and we have d{d — 1) segments to 
paint, so that 2d is actually a lower bound for the game 
turns. 

Colouring for d even - This time we draw the graph 
as a c? — 1 sided regular polygon, with the additional ver- 
tex in the geometric center, and consider all the edges, 
diagonals, and vertex-centre segments, see figure [2] At 
the first turn the player paints in red one edge, the par- 
allel disgonals, and the segments joining the center with 
the farthest vertex. Afterwards, he rotates the pattern; 
and subsequently he does the same withthe blue brush. 
When he is finished, he uses a turn to paint all the ver- 
tices. That completes the game in 2d— 1 turns. As before, 
2(i — 2 is a lower bound for painting the segments, but 
if that is the case we need at least one additional turn 
to paint the dots. Therefore, this colouring strategy is 
optimal. 



A. Multiple bounds from data processing 

In the previous section we presented a strategy for eval- 
uating (F) and (G) according to the local decompositions 



( 23 1 . Now we want to show that it is possible to repro- 



cess the same (classical) data to obtain expectation val- 
ues of other QEWs. These are generated from F and G 
by means of local invertible transformations as discussed 
in section [n] Specifically we focus on the following phase 
shift operators 



(24) 



where - ,c?}^-{ — l,l}isa generic 2- valued func- 

tion. Now, let us define 



F(0 = (1 «) W{£,)) F (1 ® W\0) 
G(0 - (1 «) W{0) G (1 (g) W^^(O) 



(25) 



For all choices of i^, these observables are QEWs of 
Werner and Isotropic class respectively (they are just 
associated with Werner and Isotropic states defined on 
a different basis). Furthermore since they are LOCC 



equivalent to the original F and G, their optimal bound- 
ing function are still given by Q. By optimizing with 
respect to all W{£,) we can write 



£of{p) > max{6wcr(Tr [F(0 p])} 
£of{p)> max{6i,o(Tr[G(0p])} 



(26) 



where the maximum is taken over all possible choices of ^. 
These expressions can be further simplified by exploiting 
the monotonicity properties of 6wcr(2;) and &iso(a;) (the 
former is non increasing while the latter is non decreas- 
ing) to move the optimization over ^ on their arguments. 
This yields 

£of{p) > bwci-if*) , £of{p) > biso{g*) , (27) 



where we defined the quantities 

Mp) 



mm{Tr[F(0 p]} 



5,(p)=max{Tr[G(e)p]} 



(28) 



These are complicated (non linear) function of the state 
p which are related to the quantity of (|6|. The main 
difference is that the latter was optimized with respect 
to the whole set of unitary LOCC transformations F, 
while and 5* are optimized only with respect to the 
a proper subset of such set, namely the phase flips (25 1. 
Luckily enough we can simplify and by rewriting 
the operators F(^) and G(^) in terms of the operators 
( 22 1 . Indeed by doing so we find that they are made of 
the very same d^ operators appearing in ( 23 ) , namely. 



F(6 = E^"®^" + iE^(")^(/^) 



X (A„/3 (g X^f3 + Y^p (^Y^p) , (29) 



and 



G(e) = E^"^^" + iE^(")^('3) 



X {X^fs ® X^a - Y^p ®Y^p) . (30) 

This implies in particular that their expectation values 
can be recovered by the same set of LOCC measurements 
required for (F) and (G), allowing us to build 2^*"^ inde- 
pendent new Werner QEWs of the form F(^), as well as 
2^-1 independent new Isotropic QEWs of the form G(^) 
by simple data post-processing [the number 2'^~^ follows 
from the fact that there are 2'^ possible ^ functions, but 
F(0 ^ F(-0 and G(0 = G(-0]- Clearly this is a very 
simple and cheap way to substantially increase the in- 
formation we can gather on p. As a matter of fact the 
expressions ( 29 1 and ( 30 ) are also useful to simplify the 



optimizations of ( 28 1 . In particular define the d 
matrix A/^*^ of elements. 



M 



(±) 



'-{{X^p (g) X^p) ± {Y^p (g) Y^p)) , a^l3 



(31) 



where, as usual, (• • • ) stands for taking the expectation 
value with respect to p. These matrices are composed by 
the measurement outcomes one acquires when recovering 
the expectation values of F and G via the local decom- 
position introduced in the previous section. In particular 
we can write, 

Tr [F(o p] = Y^ ^^p) ^ ■ • 

a, 13 
a, 13 

(32) 

where ^ = (C(l), • • • i^.id))'^ is the vector formed by the 
output entries of the binary function ^. Therefore, in 
order to determine the values ( 28 1 is sufficient to respec- 



tively minimize or maximize these quantities over ^; i.e. 



5*(p) = max e'^-Af(-) 



(33) 



Unfortunately, apart from the trivial case of d = 2, pro- 
viding a general analytical solution for these expressions 
is nontrivial and dependent on M^^^ (indeed it is for- 
mally equivalent to find a ground state of a frustrated, 
classical many-body spin system). On the other hand, 
one can easily evaluate upper bounds as follows 



/*(p)>E*e'-Ei</.'i' 

a a^f3 

5*(p)<E*^t«^ + EKii' 

a 



(34) 



even though for d > S they will be in general not tight 
and will not be useful to lower bound the entanglement 
of formation of the system. Luckily enough however the 
optimizations ( 33 ) are treatable numerically, especially 



considering that for the experimental implementation we 
are interested in, d is usually of the order of magnitude 
of units (see next section). 



IV. IMPLEMENTATION: MULTI-RAIL 
ENCODINGS IN QUANTUM OPTICS 

QEW techniques for estimating lower bounds of en- 
tanglement of a quantum state are suitable, in princi- 
ple, to any experimental setting or physical apparatus. 
Yet, we are inclined to think that the proposed local de- 
composition scheme of Werner and Isotropic witnesses 
we presented in the previous section would be particu- 
larly fit to multi-rail encoding experimental settings. In 
such implementations, a state \a) of the canonical basis 
is represented by the presence of a particle/carrier in a 
specific transmission channel (e.g. spatial mode, polar- 
ization, time-bin or orbital angular momentum) indexed 
by a. Multi-rail encodings are common both in optics 



A -< ^0 




B < ^0 



Figure 3: (colour online) Example of data sampling for 
AB system in which the local degree of freedom of each 
site is composed by three independent optical channels (here 
represented as spatially separated for simplicity). The fig- 
ure represent the fist step of the detection strategy: Ap- 
plying the 50% transmittive beam splitter and measuring 
{{Na,i — NA,2)iNB,i — Nb.2)) gives us the expectation value 
of Xi^2 <8) Xi^2 (if we put a $ = 7r/2 phase shift on channel 2, 
we would have yi,2 ® ^1,2 instead); at the same time we also 
acquire {P3 ® P3) — {Na,3 Nb,3)- For later steps, we need to 
permute the channels in the picture and perform those mea- 
surements as well. 



[27H3T] and in electronics (e.g. see [32l[33]). In the fol- 
lowing however we will focus only on the former case as it 
is intrinsically more flexible. For these architectures the 



QEW detection strategy discussed in section pll| entails 
to properly mixing a pair of incoming modes via sim- 
ple beam-splitter transformations and to perform photo- 
coincidence counting measurements at the exit port of 
the device, see figure [3] Accordingly the proposed LOCC 
measurements appear to be more suited to detect entan- 
glement in these systems than (say) the optimal schemes 
of j25j which requires coherent mixing among multiple 
channels. 

Consider for instance two separated sets A and B of 
(monochromatic) spatial optical modes. Within each of 
the two subsystems a single photon is injected and al- 
lowed to flow through d possible optical paths. In one 
of experimental setup of [25] for instance, d = 8 was re- 
alized by exploiting polarization and 2 different spatial 
degree of freedom. A generic (pure) state of the system 
can then be written as 



I*, 



E 

a, 0=1 



a/3 , 



E ^c^aP) , (35) 

a, 13=1 



where the the annihilation operator of a pho- 



7 



ton flowing through channel a at port /subsystem j G 
{A, i?}, and where \afi) stands for the two photon state 
|0). We will describe local transformations on 



photonic state as a forward elastic scattering process 



7=1 



,7 Ob, 7 

7=1 



(36) 



with U and V unitaries which acts on the A and B 
sets respectively. After applying such transformation 
all we have to do is record the channel-selective photo- 
coincindence measurements 



= \{af3\U®V\^,n)\\ 



(37) 



This way we achieve simultaneously the expectation val- 
ues of the (P one dimensional projectors of the form 
«) |a/3)(a/3| U i^V. The last step is to identify 
a proper set of unitaries U and V that lead to the de- 
termination of the quantities {Pa Pa), {Xai3 ® Xap), 
and (Yap ® Yaj}) which are needed for characterizing the 
selected QEW. For instance assume we want to measure 
the expectation value of Xap ® Xap , where a and /? rep- 
resent to two distinct channels. Then U and V could be 
implemented as Hadamard gates, acting in both subsys- 
tems A and B, which couples the pair of such channels. 
This entails the implementation of a beam-splitter trans- 
formation of transmissivity 50% which coherently mixes 
the two channels. After that, we can convert coincidence 
measurements into data via 

{Xa0 ^ Xap) = {{NA,a " ^A,/3)(^B.c< " Nb^^)). (38) 

Similarly, if our goal were to obtain (Yajj ® Yap) instead, 
the unitarities U and V could be chosen as Hadamard 
transformations, this time applied after a local phase 
shift (e.g. of 7r/2 acting upon channel /3). See figure 
[3] for details. We can now translate into experimental 
operations every possible strategy for the local QEW de- 
composition of section [Till 



As a concluding remark we point out that the same 
analysis can be generalized to an electronic implementa- 
tion of multi-rail encodings [311 133] • In this case the 
ports A and B will correspond to independent leads, 
whose transverse modes play the role of the optical paths. 
Under appropriate conditions [32] the coincidence count- 
ing ( [ST] ) can be related to the zero-frequency shot noise, 
which is a measurable quantity. In |37| . for example, a 
procedure for constructing EWs was employed for dif- 
ferent solid-state systems for the case of two transverse 
modes per lead. In such systems, however, an experimen- 
tal limitation is represented by the difficulty of working 
with more than two modes. 



V. CONCLUSIONS 

We have reviewed the notion of QEWs and clarified 
that there exist a class of such operators for which the op- 
timal bounding function can be easily determined with- 
out necessarily going through the complex optimization 
procedure of ^ and ([3]). For these QEWs we then con- 
sidered an explicit local decomposition that allows one 
to acquire the needed expectation values via a limited 
number of LOCC measurements. The proposed scheme 
is not optimal in terms of number of independent setups 
one need to prepare, but has two main advantage with 
respect to other schemes [3S]. In particular we can eas- 
ily optimize with respect a large class of unitary LOCC 
providing an operational characterization of the associ- 
ated non linear quantities, see (32 1. Furthermore we have 



shown that it is suitable for experimental implementa- 
tions based on multi-rail encodings by means of simple 
two-mode transformations (beam-splitter mixings) fol- 
lowed by correlated photo-counting detections. 
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